Shell model study of the pairing correlations 
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A systematic study of the pairing correlations as a function of temperature and angular momentum 
has been performed in the sd-shell region using the spherical shell model approach. The pairing 
correlations have been derived for even-even, even-odd and odd-odd systems near N=Z and also for 
the asymmetric case of N=Z-|-4. The results indicate that the pairing content and the behavior of 
pair correlations is similar in even-even and odd-mass nuclei. For odd-odd N=Z system, angular 
momentum 1=0 state is an isospin, t=l neutron-proton paired configuration. Further, these t=l 
correlations are shown to be dramatically reduced for the asymmetric case of N=Z-|-4. The shell 
model results obtained are qualitatively explained within a simplified degenerate model. 
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I. INTRODUCTION 

It is well established that the pairing field is an impor- 
tant component of the nuclear mean-field potential. The 
interplay between the deformation driving forces and the 
pairing field determines most of the properties of a nu- 
clear system. The relevance of the pairing field for the 
nuclear many-body system was proposed in thepioneer- 
ing work of Bohr, Mottelson and Pines in 1958 [l| . There 
are two issues related to the pairing potential which still 
need to be elucidated. The first issue concerns with the 
detailed form of the pairing potential and the second to 
the approximation employed for the solution of pairing 
force. Although, it is quite evident from the analysis of 
the properties of nuclei, for instance the suppression of 
the moments of inertia of rotating nuclei and the observed 
energy gaps, thatpairing field is essential for describing 
atomic nuclei f^, 0] . But most of the properties of nuclei 
are rather insensitive to the detailed form of the poten- 
tial in the pairing or the particle-particle channel. In 
comparison, most of the nuclear properties such as com- 
pressibility, surface energy, effective masses, saturation 
and other specific properties critically depend on the ex- 
act form of the nuclear potential in the Hartree Fock or 
particle-hole channel. In addition, the potential in the 
particle-hole channel is constrained or adjusted to the 
properties of the closed shell nuclei and the pairing field 
is introduced in an ad hoc manner in most of the density 
functional theories (DFT). In most of these approaches, 
for instance Skyrme and relativistic mean field models, 
the potential is chosen to be of zero range and therefore, 
in principle, don't contain pairing correlations. In the fi- 
nite range Gogny density functional approach, the range 
of the effective pairing force is adjusted to the properties 
of G-matrix calculated from the bare nucleon-nucleon in- 
teraction. 

The second issue concerns with the method employed 
to solve the pairing interaction. The pairing field has 
been primarily studied in the Bardeen-Cooper-Schrieffer 
(BCS) or Hartree-Fock-Bogoliubov (HFB) approxima- 
tion. In this approach, pairing field depicts a sudden 



transitional behavior as a function of rotational frequency 
and temperature 0, Q . The pairing correlations are fi- 
nite up to a certain rotational frequency or temperature 
and then these correlations suddenly vanish above this 
transitional point. The empirical analysis of the experi- 
mental data, however, does not depict this sudden phase 
transition and shows a smooth transition from one phase 
to the other. The reason for this discrepancy is known 
to arise from the iieglect of the fluctuations in the mean- 
field models 1, 0, Hi \M ■ 

It is known that relative size of the fluctuations be- 
comes small for a system with a very large number of 
particles, for example a bulk superconductor, and sud- 
den transitional behavior is a hall mark of these macro- 
scopic systems. However, for systems with flnite num- 
ber of particles, for instance atomic nuclei and metallic 
clusters, the fluctuations are important and need to be 
incorporated for an accurate description of these systems 
[Til . . A powerful method to incorporate the fluctua- 
tions is the projection theory. In particular, the particle 
number projection methods are now readily available to 
incorporate the fluctuations quite accurately. However, 
these projection methods, before variation, are available 
only at zero temperature d, [H, [13, [H, [lg| in the com- 
plete HFB framework. It needs to be added that the 
particle number projection analysis has been performed 
at finite temperature with BCS ansatz [17i, ilSj. 

The projection at finite temperature is more important 
as it is known that the mean-field or BCS wavefunction 
for an even or odd particle system has admixtures from 
both even and odd neighboring particle numbers [13, [3 • 
At zero temperature, even (odd) system has admixtures 
from only even (odd) particle numbers. It has been also 
demonstrated recently, in an exactly solvable model, that 
the pairing correlations re-appear at finite temperature 
after they are quenched at zero temperature and high 
rotational frequency (8l.[20|. This surprising result, which 
completely contradicts the mean-field predictions, needs 
to be studied in more realistic models. 

The spherical shell model (SSM) is another tool to 
study pairing correlations. The advantage of the SSM 
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approach is that, first of aU, most of the interactions 
employed in SSM are adjusted to open shell nuclei and, 
therefore, contain a proper pairing force. Secondly, the 
exact diagonalization of the Hamiltonian matrix is per- 
formed and results obtained contain all the possible cor- 
relations including the pairing. There are many puzzling 
questions regarding the pairing correlations which can be 
addressed within the framework of SSM. For instance, 
difference between the pairing correlations in odd- and 
even-mass nuclei, importance of the neutron-proton pair- 
ing correlations near N=Z and the phase transition men- 
tioned above. These issues originate when comparing the 
mean-field results with the experimental data or exact so- 
lutions of toy models. It needs to mentioned that in SSM 
approach only low-lying part of the excitation spectrum 
is fitted and, therefore, the results obtained using SSM 
are inaccurate for high-excitation energy or temperature. 

The purpose of the present work is to study the pair- 
ing correlations in a realistic space using the spherical 
shell model approach. It needs to be emphasized here 
that a new shell model program has been developed by 
two of the present authors [2l| . This new program com- 
pletely works in the j-representation and is quite sim- 
ilar in structure to that of NATHAN code developed by 
the Madrid-Strasbourg group [1^. The shell model ap- 
proach provides the most accurate description of nuclear 
properties and incorporates all the possible correlations. 
However, the SSM analysis is restricted to lighter nuclei 
and its application to heavier nuclei appears impossible 
in the near future. The most recent progress in the SSM 
approach is the study of the fp-shell nuclei [1^ . 

In the present work, the shell model analysis has been 
performed in the sd-shell region. The reason for choos- 
ing sd-shell is that, first of all, the interaction is well 
established and it is known that "USD" interaction [23| 
provides an accurate description of most of the sd-shell 
nuclei. Secondly, it is required, in the present analysis, to 
calculate a few thousand states for each angular momen- 
tum to evaluate the statistical partition function. These 
calculations would become impossible in the fp-shell re- 
gion. The calculations have been performed in the middle 
of the sd-shell so that a large number of eigen-solutions 
are available for the statistical analysis. The nuclei for 
which the shell model calculations have been performed 
are : ^^Si, ^^Si and ^^Al. These three neighboring nuclei 
have been investigated so that a comparison can be made 
among even-even, even-odd and odd-odd systems. The 
pairing correlations have also been deduced for the asym- 
metric system ■^'^Ne in order to study the dependence of 
the neutron-proton pairing on the particle number. 

It is pertinent to mention here that while the present 
work was in progress, a similar study of pairing correla- 
tions in the sd-shell has been recently published [1^ . Al- 
though the model and the region of study in the present 
work is same as that of ref. [1^, the issues discussed 
are different and as a matter of fact the present work 
complements the work of ref. [l^ . In the course of the 
discussion, we shall comment on the results of ref. 1291 



wherever necessary. The present manuscript is organized 
as follows : In the next section, the shell model based ex- 
pressions are presented for the evaluation of the pairing 
correlations. The results of the calculations are presented 
and analyzed in section III and finally the summary and 
conclusions are included in section IV. 



II. SHELL MODEL FORMULATION 

The spherical shell model Hamiltonian, generally, con- 
tains single-particle and two-body parts and in the sec- 
ond quantized notation is written as 



H = hsp + V2, 



where, 
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where e^s are the single-particle energies of the spheri- 
cal shell model states, which are diagonal except in the 
radial quantum numbers and < rs\va\tu >r are the two- 
body interaction matrix elements and in the present work 
are chosen to be those of "USD". The two-particle cou- 
pled operator in Eq. (3) is given by ylp(rs) — (cjcj)r 
and ArMr = (— l)'"^^^'^^r-Mr • The labels r, s, ... in the 
above equations denote the quantum numbers of spin 
and isospin. "F" quantum number labels both angular 
momentum and isospin of the coupled state. The above 
notation is same as that used in ref. [2^ . 

In the present work, the pairing correlations have been 
calculated using the canonical ensemble approach since 
the exact solutions have well defined particle number. 
The average value of a physical quantity "F" in canonical 
ensemble is given by 0, 



«F»=J2F^e'^'/''^/Z, 



where, 
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F, = <i\F\i> 
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In the partition function, Eq. (5), k is the Bohzman 
constant and T is the temperature. In the rest of the 
manuscript, we shaU use "Temp" rather than kT, which 
has dimensions of energy (MeV ) and the symbol T shall 
be used for the isospin of the coupled two-particle state. 

In the present work, we study the isovector monopolc 
pair correlations for the 'Tq = {/ = 0,T — l,\Tz\ = 
0, 1}". The canonical average for this is calculated as 

/ N v'(2r + 1) , , 

^^ro(pair) = 2^ === < rs\va\tu >ro 

rstu V(l+'5rs)(l+dt„) 



<< {Al^ X Ar„)o », 



(7) 



from the energies and eigen-states obtained by diagonal- 
ization, which shall be referred to as the "correlated" 
pairing energy Epair{corr). We have also calculated the 
"uncorrelated" pairing energy by removing the monopole 
field Fq in the interaction and then redoing the shell 
model diagonalization Q. Using the resulting energies 
and wavefunctions, Eq. (7) is evaluated to give the "un- 
correlated" pairing energy. It is noted that the definition 
of present pairing energy is slightly different from that 
used in refs. [25l 12^. In these studies, it is calculated 
without the matrix elements and the constant factors in 
front of the expectation values in Eq. (7). 



III. RESULTS AND DISCUSSIONS 

The shell model calculations have been performed in 
the middle of the sd-shell for ^^Si (with six valence pro- 
tons and six valence neutrons), ■^^Si (with six valence 
protons and five valence neutrons) and ^^Al (with five 
valence protons and five valence neutrons). The calcula- 
tions have also been performed for ^''Ne (with two pro- 
tons and six neutrons) in order investigate the behav- 
ior of the pairing correlations with asymmetry in proton 
and neutron particle numbers. The temperature and the 
angular momentum dependence of the pair correlations 
have been studied in detail for these systems and are dis- 
cussed in subsections B and C. In subsection D, the tem- 
perature dependence of the average isospin is discussed. 
Before presenting these results, in subsection A a simple 
model is briefly described to analyze the essential results 
obtained with the full shell model calculations. 



A. Isospin geometry and a simple model 

Before starting the discussion, we first briefly recall 
some consequences of isospin being a good quantum num- 
ber. Writing explicitly the isospin and omitting the other 
quantum numbers, i.e. T — T,Tz, the pair energy is mea- 
sured by 
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xii,_M= <^Ml- M\tO> Bto, (8) 

t=0,l,2 



where the proton-neutron paring corresponds to M = 
and the neutron- neutron pairing to M = I. Since A J 

and Ai -M are tensor operators in isospace their product 
can be rewritten as a sum of tensor operators Bto and 

< TT,\Al^,j X ii,_M|Tr, > 

= ^ < IMl - M\tO > 

t=0,l,2 

X < TT;,10\TTz >< T\\Bt\\T > . (9) 

Only the t = term contributes in T = states, which 
means that < T = 0|AJ x Ai^^m\T = > does not 
depend on M, i. e. Epair is the same for all three types of 
pairing. For T > states, the terms i = 1, 2 contribute as 
well, which means Epair may be different for each pairing 
channel. 

For a quantitative statement, one needs to know the 
reduced matrix elements < T||i?(||T >, which depend on 
the detailed structure of the shell model states. For the 
lowest states, pairing correlations are strong. In such a 
case one expects that the model of isovector monopole 
pairing in a degenerate shell should allow us some rough 
estimate of the relative paring strengths. The model is 
discussed in [13, [13 1 where the original work is cited 
and explicit expressions are given in terms of number of 
nucleon pairs Af in the shell, the isospin T,Tz, and the 
number of pairs ri(= 12) that the shell can accommodate, 
where the unit is the coupling constant G. It was found 
there that the model accounts well for relative strengths 
of the pair energies but cannot reproduce the scale of 
the shell model calculations, which was attributed to the 
splitting of the levels due to the deformation and the 
spin-orbit coupling. Following this observation, we scale 
the pair energies of the degenerate model by a common 
factor. 



B. Temperature dependence of the pair 
correlations 

The results of the neutron- neutron and neutron-proton 
monopole pair correlation energies for f 48114 are shown 
as a function of temperature (Temp) in Figs. 1 and 2 
for even- and odd-spin values. For the symmetric system 
and isospin invariant two-body interaction, the proton- 
proton pairing energy is identical to that of neutron- 
neutron. The reason is that the canonical ensemble con- 
tains mainly isospin, T = states, which lie lower than 
the T > states in the N=Z=even nuclei (29| . 

Even- and odd-spin values are plotted separately in 
two figures since they have different intrinsic structures. 
In the quasi-particle language, the low-lying even-spin 
members in an even-even nucleus have 0-quasiparticle in- 
trinsic structure and the odd-spin members have 2-quasi- 
particle structure. Obviously, the fraction of T > states 
is larger for odd spin. 

In order to investigate, in detail, the variation of the 
pairing correlations with spin (/), the correlations are 
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plotted separately for each possible value of spin and 
are shown in two panels. In the top panel the pairing 
correlations are calculated with full two-body interac- 
tion [denoted by Epair{corr)]. The lower panel depicts 
pairing correlations [denoted by Epairiuncorr)] in which 
the monopole terms in the two-body interaction were ex- 
cluded. 

For / = in Fig. 1 , the Epair (corr) are quite similar for 
neutron-proton and identical particle channels and are 
also quite large at low temperatures. For Temp=0, they 
are 4.6 MeV. The ground-state in even-even systems is 
a paired configuration with maximum correlations. The 
degenerate model gives 3/2 for the correlation energy, 
which fixes the scale factor to 3.06 MeV. In the following 
we shall use this factor to scale other calculations in the 
framework of the degenerate model. 

The pair correlations are almost unchanged till 
Tenip~2 MeV and then are observed to be reduced with 
increasing temperature. However, in comparison to the 
mean-field models which predict a sudden transition from 
the paired to the unpaired state, the exact analysis de- 
picts a smooth drop in the pair correlations. The phase 
transition obtained at Temp~2 MeV in Fig. 1 is higher 
as compared to the transition point obtained in HFB 
and also in shell model Monte-Carlo study, which pre- 
dict at Temp~0.5 to 0.7 MeV It is to be noted 
that in the present work, the pairing correlations have 
been obtained using the expression, Eq. (7) and, as al- 
ready pointed out at the end of section II, is slightly 
different from the expression used in the HFB and the 
shell model Monte-Carlo studies. In order to confirm 
that the reason for the discrepancy in the phase transi- 
tion point is due to the different pairing expressions used, 
we have performed shell model calculations for a simpler 
case of ^^Mg with the pairing expression as used in HFB 
and Monte-Carlo studies and the phase transition was 
found to be around 0.9 MeV. It is also noted that the 
magnitude of the pairing energy obtained in the present 
study is also quite large as compared to the HFB and 
the shell model Monte-Carlo results. Further, it is seen 
in Fig. 1 that at large temperatures the neutron-proton 
pairing energy deviates from proton and neutron pairing 
energies, which was also observed in ref. [2^. The dif- 
ference between the proton-neutron and neutron-neutron 
(=proton-proton) pair energies indicates some admixture 
of r = 1 states to the ensemble. 

For even-spin values of 1—2, 4 and 6, it is observed that 
the pair correlations drop in a step wise manner and the 
pairing gaps for these spin values at Temp=0 are approx- 
imately 3.9, 3.2 and 2.6. The temperature dependence of 
the pair energies for 1 = 2 and 4 show a similar behavior 
as that of / = 0. For 1=6, the pair correlations are almost 
constant with temperature. The shell model calculations 
for 7 = 8 and 10 (not shown in the figure) depict a slight 
increase in the pair correlations at low temperatures and 
for higher temperatures the correlations drop as for the 
earlier cases. 

The uncorrelated pairing energy are shown in the lower 



panels of Fig. 1. As already mentioned, they have been 
obtained by a second shell model diagonalization on set- 
ting the monopole matrix elements equal to zero. The 
problem in the shell model study of the pair correlations 
is that the calculated correlated pairing energy may con- 
tain the contributions from other multipoles due to re- 
coupling, which in the mean-field language are referred 
to as the particle-hole contribution. In the mean-field 
framework, particle-particle and the particle-hole chan- 
nels are completely decoupled and the pairing energy can 
be directly evaluated from the pairing potential. 

It is noticed from Fig. 1 that the uncorrelated pair- 
ing energy is substantially smaller as compared to the 
correlated one. However, it shows a similar transitional 
behavior with temperature as the correlated pairing en- 
ergy. In order to explore the reason for this unexpected 
transitional behavior in the uncorrelated energy, we have 
also removed 1=1 and 2 apart from 1=0 matrix elements 
in the shell model diagonalization and the pairing energy 
again calculated using the expression of Eq. (7). The 
results indicate that the phase transition is now almost 
washed out in the uncorrelated pairing energy. We are 
presently performing a detailed investigation of this phe- 
nomenon and the results of this study will be presented 
elsewhere. 

The neutron-proton uncorrelated pairing energy ap- 
pears to be lower than that of the corresponding neutron 
and proton energies. This can be understood as follows. 
The isovector monopole pair correlations shift the T = 
states to lower energy as compared to the T > states. 
If these correlations are switched off, then this prefer- 
ence is lifted. The increased fraction of T > states in 
the ensembles creates the difference between the neutron- 
neutron and proton-neutron pairing. However, using this 
simple perspective, it is not possible to understand why 
the like-particle pairing is stronger. 

The odd-spin values are depicted in Fig. 2. The pair 
energies are smaller than the even-spin values at very 
low temperature. This is easily understood by noting 
that the odd-spin band have two quasi-particle structure 
with reduced pairing. For / = 1, the yrast state has T = 
0. The pronounced increase of the neutron-proton pair 
energy as compared to the like-particle ones indicates 
that low-lying T > states become a substantial fraction 
of the ensemble. 

For the larger angular momenta, / = 7 it is observed 
that the pair correlations first increase with temperature. 
It is well known that in the mean-field theory, the pair 
correlations are reduced with both increasing temper- 
ature and rotational frequency. The mean-field analy- 
sis always depicts a phase transition from the paired to 
the unpaired phase with increasing temperature and ro- 
tational frequency. This is experimentally observed in 
macroscopic superconductors, where the particle num- 
ber is quite large. However, for mesoscopic systems the 
mean-field solution is inappropriate. The exact solutions 
of the pairing problem do not have a sharp phase tran- 
sitions. Further, they may have peculiarities, as e.g., 



5 



recently demonstrated in a simple few-level model that 
the pairing correlations re-appear at finite temperature 
after they are quenched at zero temperature and high 
rotational frequency where this peculiar effect is ex- 
plained. The shell model results for / =9 and 11 (not 
shown in the figure) also depict a minor increase. There- 
fore, increase in the pairing correlations with increasing 
temperature obtained in the earlier work is an artifact 
of the simple model employed as the present results 
in the sd-shell don't depict such a dramatic rise. How- 
ever, it needs to be added that in the sd-shell, it is not 
possible to have a large angular momentum where the 
re-appearance of the pairing correlations is predicted in 
the earlier simple model study. 

The results of the pair energies for the odd-system 
ilSiia are presented in Fig. 3. The results presented 
in this figure also correspond to 13AI14 with proton and 
neutron curves interchanged. At zero temperature the 
expected picture evolves. The proton-proton energy of 
4.4 MeV is similar to the value of 4.6 MeV in the even- 
even neighbor ^48114, as both have same proton number. 
The neutron-neutron pair energy of 3.0 MeV is about 
30% lower due to the blocking of one level in the odd- 
neutron system. This is in accordance with the mean- 
field predictions that the pair gaps in an odd-system 
is reduced by about 15% in comparison to the neigh- 
boring even-even system, which amounts to a reduc- 
tion of about 30% in the correlation energy. The np- 
pair energy of 3.7 MeV is less reduced than the neutron- 
neutron one. This is understood, because the odd neu- 
tron blocks the level only partially for proton-neutron 
pairs. [If the odd neutron is in the state m > 0, a pair 
with the neutron in m < and the proton in to > 
may by accommodated.] The pair energies are differ- 
ent, which is expected for T < 1/2 states from Eq. 
The degenerate model allows us to estimate the conse- 
quences of blocking. For T = 1/2, TV = 4 it gives 
Epp{corr) = 3/2, E„p{corr) = 15/12, Enn{corr) = 1. 
The respective scaled values of 4.6, 3.8 and 3.1 MeV com- 
pare well with the shell model values of 4.4, 3.7 and 3.0 
MeV. 

In early nineties, it was demonstrated through a series 
of systematic studies [3l| of experimental data that the 
difference between the moments of inertia of neighboring 
even-even and odd-nuclei is merely < 2%. This posed 
a serious challenge to the traditional nuclear structure 
models based on mean-field theory which give a differ- 
ence of about 15%. It was shown later using number 
projected mean-field models that this difference in the 
moments of inertia could be reduced. The moment of 
inertia depends on deformation and the pairing corre- 
lations and considering that ^^Si and ^^Si have similar 
deformation values, it is expected that the moments of 
inertia of two nuclei would be similar since the pair gaps 
of the two systems are similar. The calculated moments of 
inertia for low-lying states for ^*Si and ^^Si are 2.62205 
fi^/iMeV) and 2.66220 fi^/iMeV), respectively. 

The results for the odd-odd fsAlis system are pre- 



sented in Figs. 4. For low-spin, the pair energies are 
rather different from those of the even-even and odd- 
systems. For / = 0, neutron-proton pairing energy is 
quite large as compared to like-particle pairing energies. 
The ground states in odd-odd self-conjugate nuclei show 
a preference for T = 1 [2^ with a T = state close by. In 
our case the ground state has T — 1. For T — 1, — 0, 
the simple model of isovector pairing in a degenerate 
shell gives a ratio Epn{corr) — 43/15, Enn{corr) = 
Epp{corr) = 11/15. The respective scaled values of 7.8 
and 2.2 MeV are to be compared with the shell model 
results of 7.9 and 2.7 MeV. The proton-neutron pair en- 
ergy remains larger than the neutron-neutron pair energy 
for finite temperatures, which indicates that the isovector 
proton- neutron pair correlations lower the T = 1, = 
states relative to the T = 0, = states. The lowest 
1 = 2 state has T = and the same pairing energies 
for like and unlike particles. For higher temperature, 
neutron-proton pairing is larger as compared to identi- 
cal particle pairing, which indicates a large fraction of 
T = 1 states in the ensemble. With increasing spin the 
results become similar to the even-even system, which 
can be understood as a consequence of the quenching of 
the pair correlations in both types of nuclei. This is cor- 
roborated by the observation for the odd-A case that the 
different pair energies become similar (though not equal) 
at large spin. For odd-spin values, the pair energies are 
nearly equal for all spin values. Only at higher temper- 
atures there is a slight asymmetry. This indicates that 
the odd-spin states have preferentially T = 0. 

The pair correlations for the asymmetric system ^QNei4 
are presented in Fig. 5. As expected, the neutron pair 
energy is larger than the proton one, because there are 
more neutrons in the open shell. The neutron-proton pair 
energy is quite small. In contrast to the symmetric even- 
even system, where it is as large as the like particle one. 
This is in accordance with the HFB theory with neutron- 
proton pairing 32] , which finds quite strong neutron- 
proton pairing for the N=Z system, but vanishing for 
the asymmetric case of N=Z-|-4. The present shell model 
substantiate these HFB results, although the neutron- 
proton has not vanished in the shell model study but 
has clearly become quite weak. The reason is that the 
extra neutron pair blocks a level for the proton-neutron 
pairs to scatter into. For zero temperature, the difference 
between the symmetric and asymmetric systems is qual- 
itatively reproduced by the simple model of a degenerate 
shell. The ground state has A/" = 4, T = 2, which gives 
Epp{corr) = 13/14, E„n{corr) = 27/14, Enp{corr) = 
13/42. The respective scaled values are 2.8, 5.9 and 1.0 
MeV to be compared with the shell model results of 3.4, 
5.5 and 1.2 MeV. 

The model of a degenerate shell obeys the dynami- 
cal symmetry of the group SO (5) with the isospin sub- 
group SU(2). As will be discussed in a forthcoming pa- 
per, the relative strength of the three types of isovector 
pairing mainly reflects the geometry of the isospin in- 
duced by SU(2), which remains valid for the case of a 
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non-degenerate shell. The full SO (5) symmetry is ex- 
ploited when the effects of blocking in states with non- 
zero seniority are estimated as for the case of odd-A in 
this paper. 



C. Angular momentum dependence of the pair 
correlations 

The angular momentum dependence of the pair corre- 
lations for the four systems studied are depicted in Fig. 6. 
The results for ^*Si indicate that for zero temperature, 
the pair correlations drop monotonically with increas- 
ing angular momentum, except for / = 1 which shows a 
larger drop. In the case of ^^Si, the pair correlations show 
a staggering effect, where the phase of the staggering is 
same for the three pairing modes. 

For ^^Al, neutron- proton correlations are maximal for 
/ = 0. This is evident from Fig. 6 with a drop of about 
5 MeV from / = to 1 in the correlated pairing energy. 
For the higher spin values, it is noted that pair corre- 
lations drop very little. The results for ^''Ne depict a 
larger staggering effect as compared to ^*Si for identical 
particle pairing. The neutron-proton pairing shows an 
irregular behavior with spin for this asymmetric system. 
As expected, for higher temperature of 3 MeV, the pair 
correlations depict a smoother behavior with spin. 



D. Isospin analysis 

In most of the analysis presented in the above subsec- 
tions, the isospin content of the states played a crucial 
role to understand the behavior of the pairing energies. 
It is, therefore, imperative to ascertain the temperature 
dependence of the isospin. It needs to be mentioned 
that the new shell model program developed [2l| uses 
the neutron-proton product basis and, therefore, isospin, 
although conserved, needs to be evaluated for each eigen- 
state. The expectation value of has been calculated 
using the shell model wavefunctions as discussed in ref. 
[3^ 1 . The average value of as a function of tempera- 
ture has been obtained by using the canonical partition 
function in the same manner as the pairing correlations 
have been deduced. 

The average value of is shown in Fig. 7 for the 
lowest angular momentum ensemble. In the top panel 
of Fig. 7, average value of is plotted for ^^Si and it 
is quite evident from this figure that the isospin at low 
temperature is T=0 up to Temp=2 MeV. Above this tem- 
perature, the isospin increases steadily and correlates well 
with the drop in the pair energies observed in Fig. 1. For 
the odd-odd ^^Al system, the isospin at low temperature 
is equal to one and giving rise to large difference in the 
pairing correlations between neutron-proton and identi- 
cal particle channels in Fig. 4. However, with increasing 
temperature, it is noted that the average isospin drops as 
T=0 states enter into the ensemble with the consequence 



that neutron-proton pairing energy comes closer to the 
identical particle pairing at higher temperatures in Fig. 
4. 

For the asymmetric case of ^"^Ne, the average isospin 
in Fig. 7 is equal to two and appears to be almost con- 
stant with increasing temperature. This constancy of the 
average isospin gives rise to the constant behavior of the 
neutron-proton pairing energy in Fig. 5. The tempera- 
ture behavior of the isospin for the odd-mass system ^^Si 
is similar to the even-even system ^^Si and is the reason 
that the behavior of the pairing correlations for the two 
systems are quite similar. 



IV. SUMMARY AND CONCLUSIONS 

In the present work, the shell model study of the pair- 
ing correlations has been undertaken. The calculations 
have been performed in the sd-shell for ^^Si, ^^Si, ^^Al, 
and ^"^Ne. For the case of the even-even system ^^Si, the 
pair energy of even-spin states as a function of temper- 
ature depicts a smooth but pronounced decrease around 
Temp=2 MeV, which can be interpreted as the strongly 
washed out relic of the phase transition. In the case of 
the odd-spin states up to / = 7, the pairing correlations 
decrease only very slowly with increasing temperature, 
starting from a reduced value, which is caused by the 
blocking of two levels in the states carrying two quasi 
particle character. 

It is also clearly evident from the present study that 
pairing correlations are non-zero even at large tempera- 
tures and angular momenta. This is in contradiction to 
the mean-field predictions that pairing correlations die 
out at higher temperatures and angular momenta. For 
the very small systems studied, most of the pair corre- 
lations is generated by the fluctuations of the pair field, 
which are more prominent than the mean-field itself. 

The proton-proton, neutron-neutron, and proton- 
neutron isovector pair correlation energies are not scalar 
under rotation in isospace, which means that they are 
only equal in T = states but generally different in 
T > states. For the ground states of evcn-A nu- 
clei, the relative strengths of the different pairing en- 
ergies is qualitatively reproduced by the simple model 
of a degenerate shell. The even-even N — Z nu- 
cleus has a T =0 ground state and thus equal proton- 
proton, neutron-neutron, and proton-neutron strength. 
The proton-neutron strength deviates from the identi- 
cal ones with increasing temperature because T > 
states enter the ensemble. The odd-odd N = Z nu- 
cleus 13AI13 has a r = 1, Tz — ground state for which 
the proton-neutron strength is about three times larger 
than the like particle strength. With increasing temper- 
ature, T = states enter the ensemble, reducing the the 
proton-neutron contribution. The asymmetric nucleus 
ioNei4 has a T = 2, = 2 ground state with the proton- 
neutron pairing dramatically reduced as compared to the 
symmetric systems ^^Si and ^^Al. In the odd- neutron 
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nucleus i^Siis, the proton-proton strength is about the 
same as in ^48114, the neutron- neutron strength is about 
30% lower and the proton-neutron strength in between, 
which reflects the blocking of one level by the odd neu- 
tron. 

Finally, we would like to mention that the results of 
the present work are questionable at higher temperatures 
as the configuration space of sd-shell employed in the 
present work is suited only for low-excitation energies. 
For higher temperatures, it is expected that fp-shell will 
be populated and for accurate evaluation of pairing cor- 



relations it is essential to include fp-shell configuration 
space in the shell model analysis. However, it is impossi- 
ble to perform shell model calculations with a complete 
sdfp-configuration space. What is feasible is to calcu- 
late the partition function of the fp-shell in a spherical 
degenerate limit and then calculate the total partition 
function with the method illustrated in ref. ^3^]. We are 
presently working to evaluate the pairing correlations us- 
ing this approach and the results of this analysis shall be 
presented in the near future. 
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FIG. 1: (color online) Temperature dependence of pair correlations for Si. The results are shown for even-spin values of I 
= 0, 2, 4, 6. For each spin state, there are two panels. The upper panel depicts the pairing correlations calculated with the 
full two-body interaction. The lower panel shows the pairing energy without monopole matrix elements, referred to as the 
uncorrelated contribution. For 1=0, the temperatures of 1,2, 3, 4 and 5 MeV correspond to excitation energies of 0.04, 0.90, 
4.0, 9.57 and 15.38 MeV, respectively. 
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FIG. 3: (color online) Temperature dependence of pairing energy for ^■'Si. The results are shown for I = 1/2, 5/2, 9/2 and 13/2 
on the upper part and for I = 3/2, 7/2, 11/2 and 15/2 on bottom panels. 
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FIG. 4: (color online) Temperature dependence of pairing energy for ^''Al. The pair-gaps are plotted for I = 0, 2, 4, and 6 on 
upper two panels and for 1=1,3, 5, and 7 on the lower two panels. 
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FIG. 5: (color online)Temperature dependence of pairing energy for 2*Ne.The pair-gaps are plotted for I = 0, 2, 4, and 6 on 
upper two panels and for 1 = 1,3, 5, and 7 on the lower two panels. 



